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Abstract. 

Semiclassical multimode laser theory is extended to gas lasers with open two- 
dimensional resonators of arbitrary shape. The Dopplcr frequency shift of the linear- 
gain coefficient leads to an additional linear coupling between the modes, which, 
however, is shown to be negligible. The nonlinear laser equations simplify in the special 
case of wave-chaotic resonators. In the single-mode regime, the intensity of a chaotic 
laser, as a function of the mode frequency, displays a local minimum at the frequency 
of the atomic transition. The width of the minimum scales with the inhomogeneous 
linewidth, in contrast to the Lamb dip in uniaxial resonators whose width is given by 
the homogeneous linewidth. 



PACS numbers: 42.55. Zz, 42.55.Ah 
1. Introduction 

Theory of gas lasers was developed some decades ago and is well covered in the 
literature [1,2]. The gas lasers are characterized by their capability to emit light at 
frequencies that differ from the atomic-transition frequency by more than a homogeneous 
linewidth. This property is a consequence of the coherent amplification by moving 
atoms that interact with the electric field at Doppler-shifted frequencies. If a resonance 
frequency of the cavity is continuously changed, the intensity of the single mode drops 
within the homogeneous linewidth of the atomic frequency. This narrow local minimum 
in the inhomogeneously broadened line, the Lamb dip [3], can be observed in uniaxial 
cavities that support standing waves. The Lamb dip has been utilized for stabilization 
of laser frequency [4] and precise detection of atomic-transition frequency and line shape 
(Lamb-dip spectroscopy) [5]. 

The theory of lasers (including gas lasers) was originally intended for low-loss 
resonators with a well-defined direction of propagation. The continuing interest in 
random [6, 7] and chaotic [8] lasers has motivated extensions of the conventional theory 
to irregular and strongly open systems [8-13]. These improvements relate to a number of 
aspects, such as construction of a quasimode basis of the open cavity, description of linear 
and nonlinear coupling between the lasing modes in irregular lasers, and incorporating 
the statistics of eigenfrequencies and eigenfunctions in the laser equations (see Ref. [14] 
for a review). 
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In the present article, the semiclassical theory of gas lasers is generalized to lasers 
with open two-dimensional resonators of arbitrary shape. Equations for the amplitudes 
and frequencies of lasing modes are derived by expanding the atomic polarization in 
the powers of electric field and keeping the linear and third-order contributions. Since 
the dependence of the linear gain on the frequency is modified by the Doppler shifts 
in the frames of moving atoms, there exists an additional linear coupling between the 
modes. This effect is similar to the gain-induced coupling in systems with spatially 
inhomogeneous refractive index [11] or pumping [14]. We can argue, however, that the 
the magnitude of the coupling in gas lasers is negligible under reasonable assumptions. 
The nonlinear terms in the laser equations simplify for the weakly open cavities of 
chaotic geometry, where the eigenmodes can be modeled by a superposition of plane 
waves propagating in random directions [15, 16]. In this system, the intensity of a single 
lasing mode, as a function of its frequency, has a local minimum at the frequency of the 
atomic transition. In contrast to the Lamb dip in uniaxial resonators, the width of the 
minimum is comparable to the inhomogeneous linewidth. 



2. Semiclassical theory of gas lasers 

2.1. Semiclassical laser equations 

We consider an open resonator of arbitrary shape, which, in general, can be filled 
with a homogeneous inactive medium characterized by a real dielectric constant e. 
The lasing action is provided by active atoms located inside the resonator. In the 
simplest model, only the the two atomic levels responsible for the lasing transition are 
described dynamically, whereas the information about the remaining degrees of freedom 
is contained in the empirical parameters: the homogeneous atomic linewidth 7j_, the 
population- inversion decay rate 7y , and the unsaturated population difference (per unit 
volume) An . Then the system of atoms can be characterized by the polarization P(r, t) 
and the population-inversion density (the difference between the populations of the 
upper and lower levels per volume) An(r,t). The interaction of these quantities and 
the classical electric field E(r,t) in the resonator is specified by the semiclassical laser 
equations [2] 

e d 2 E _ _ 4tt 9 



+Vx(VxE) = -^P M) (1) 
S + 2i£ + v 2 )P = -2^E(r, t) An(r, t), (2) 



dt 2 '^dt J h 

<9An , . . , 2 . dP(r,t) 

_- 7 ,(A nc -A n ) = -B(r, t )~i i li, (3) 

where v is the atomic-transition frequency and d is the atomic dipole matrix element. 
The first equation follows from the Maxwell equations (written in the Gaussian units) 
and has the form of a wave equation with the polarization as a source. The second 
equation describes the dipole interaction with the field and the decay of polarization 
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due to stochastic degrees of freedom. The third equation describes the energy exchange 
between the atoms and the field, as well as the relaxation of the population inversion. 
If the field were absent, An(r,t) would reach the stationary value of Ano- We will use 
this parameter as a measure of the pump strength. 

In the gas laser, the active atoms are moving with constant velocities Vj, which 
are assumed to be Maxwell distributed. The polarization and the population-inversion 
density can be expressed in terms of the dipole moments Pj(t) and the population 
differences ANj(t) of the individual atoms as 

M 

PM) = $>,(tH(r-r,(t)), (4) 

3=1 

N 

An( r ,t) = J2^N J (t)5(r-r J (t)), (5) 

i=i 

where M is the number of atoms and rj(t) = r° + Vjt is the trajectory of the jth 
atom. The initial positions r° are uniformly distributed within the resonator. The 
change of direction due to collisions or reflections at the boundary can be neglected if 
the time between these events is longer than the longest time scale in Eqs. (l)-(3), i.e., 
l,L ^> v/j\\, where / is the mean free path, L is the size of the system, and v is the root 
mean square velocity. 

We will consider the case of a quasi-two-dimensional resonator, and assume that 
the field vector E(r, t) has a fixed polarization direction perpendicular to the plane 
of the resonator. This assumption is consistent with the Coulomb gauge condition 
V • [eE(r,t)] = 0. Then the electric field and the dipole moments can be described in 
terms of their out-of-plane components, E(r,t) and Pj(t). 

Clearly, Eqs. (2) and (3) are local in the sense that the position r enters only 
as a parameter. Hence, the contributions of individual atoms can be separated by 
substituting the sums (4) and (5) and collecting the terms proportional to the respective 
delta functions: 

^ + 2 1± - + v 2 ) Pj = -2u-E(r 3 (t),t) AN 3 (t), (6) 
dAN j /\ at . a_ 7 \ 2 , . . , dp-jit) 

where AN = An V/M is the unsaturated inversion per atom and V is the resonator 
volume. 

Since the field in Eqs. (6) and (7) is evaluated at the positions of moving atoms, it 
is convenient to perform the Fourier transform 

£M) = 7^2 / dqE^e^, (8) 



(2tt) 2 j 

E n (t) = J^drE(r,t)e-^\ (9) 
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where the vectors r and q are two-dimensional. The transformation is applied only to the 
field inside the cavity, the integration being performed over the cavity interior X. In the 
moving frame of the atom, the Fourier component Eq(t) is multiplied by an additional 
oscillating factor exp[iq • Vjt], which expresses the Doppler shift for the frequency. In 
order to facilitate subsequent handling of lasing field oscillations, we carry out the 
temporal Fourier transform as well: 

Eq(t) | I /.',,„ J 

Pi(t) \ = 7 r dul p JW \e-**, (10) 
AA^(0 J ~ nJ - I AN* J 

(11) 
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In the (q, oj) representation, the atom equations (6) and (7) assume the form: 
(-u 2 - 2i^x_u + v 2 ) p jbJ 

= - 2 »J-J^3 J dqdw'e^E^ANj^,^., (12) 
'-iu + ^ANju = 27t 1ii AN 5(lo) 

J dqduj'u'e 1 ^ pj^E^-uji+n-vy (13) 



2i 1 



hu (2tt) 3 

Together with the wave equation (1), the above equations form a closed system, which, 
in principle, allows one to determine the electric field in the cavity. 



2.2. Third- order perturbation expansion 

In order to derive an equation for the electric field alone, the dipole moments 

Pjui can 

be expressed in terms of the field using Eqs. (12) and (13) and then substituted in the 
polarization term of the wave equation (1). To eliminate the inversion variable ANj^ 
from the atom equations, we expand pj^ and ANj^ in powers of the field amplitudes E nuJ . 
Here, the perturbation expansion is carried out up to the third-order terms, which 
provides the minimal description of the nonlinear effects and is a standard procedure in 
the laser theory [1, 2]. The third-order theory becomes invalid if the mode intensity is far 
above the lasing threshold, and a fully nonlinear approach along the lines of Refs. [12, 13] 
has to be developed. The nonperturbative theories rely on the possibility to group 
terms according to their order of magnitude in the small parameter 711 /Sou, where 5u 
is the frequency spacing between the lasing modes. These contributions can be traced 
back to oscillations of the population inversion, the zeroth order corresponding to the 
stationary- inversion approximation of Ref. [12]. The oscillations of the population at the 
beat frequencies in the multimode regime yield a correction to the polarization (linear in 
J\\/Suj and third order in the field), which oscillates at the same lasing frequency as the 
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main contribution (zeroth order in ^/Soj) to the polarization [1, 17]. This is the level 
of approximation used in the present section. The leading-order corrections in 711 /5a; to 
the fully nonlinear theory of Ref. [12] were derived and numerically studied in Ref. [18]. 
A diagrammatic expansion in 711 /5a;, nonperturbative in the field, was proposed in 
Ref. [13]. In the frequency representation, the inversion pulsations are described by the 
function D\\{ui) defined below. In the laser with gaseous active medium, the Doppler 
frequency shifts mix the orders of magnitude of the pulsation terms, and additional 
procedures have to be introduced in order to approximately reclassify them. In Sec. 3.2, 
the stationary- inversion approximation is obtained for the nonlinear contribution of the 
third order in the electric field. 

We start with the zeroth-order approximation to the population inversion by setting 
it equal to the unsaturated value: AN^j = 27rAiVo 5(o»). Substitution in Eq. (12) yields 
the linear approximation to the dipole moment, 

P# = -i^AN D(u)j±y 2 J dcle^ £ qV+q , Vj , (14) 

where we introduced the notation 

v l — uj 1 — 2i^x_uj 

The real part of D(oj) provides the frequency dependence of the linear gain for atoms at 
rest. By substituting pf^ in Eq. (13), we obtain the quadratic correction to the inversion, 



(2) 

ANjJ, which is substituted back into Eq. (12) to produce the nonlinear contribution, 



j 

72 \ 2 



P$ =jj^ (^) AiVo D(co)j±y 8 J d^d^'d^" ( hj'duj"u'D(u') D {l (u - J' + q'" • Vj) 
xe *(q+q +q ) r j^ q , jW _ w ,_ w ,, + ( q , +q ,,,). Vj E q // ia ,/ +q //. Vj E q /// w //, (16) 

where 

The resulting expansion for the dipole moment, pj^ = pfj, + pf^, serves as a 
source in Eq. (1), which, after taking into account Eq. (4) and transformation to the 
u representation, takes the form 

(eo; 2 + c 2 V 2 )^(r) 

= (2^2/ dqe^-^iuj-a. Vj ) 2 Pj> _ q . vr (18) 

A possible means to solve this differential equation is to search for the solution £^( r ) 
in the form of an expansion in the eigenfunctions of the differential operator. In the 
case of an open cavity, one can either employ the continuous basis of eigenmodes of 
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the infinite system or formulate a non-Hermitian eigenvalue problem by separating 
the functional spaces of the interior and exterior of the cavity. In the latter case, 
the complex eigenfrequencies provide the oscillation frequencies and lifetimes of the 
decaying quasimodes, whereas the right eigenfunctions yield their field distributions. 
Here we make use of the quasimode basis, as the quasimodes supply, in most cases, the 
correct description of lasing modes near the lasing threshold. 

Specifically, we employ the constant-flux (CF) quasimodes [12], which are 
particularly suitable for laser problems (see Appendix A for more information). The 
CF modes are characterized by smooth wavefunctions ipku( r ) depending on the real 
continuous parameter uj. In the asymptotic region outside the cavity, ipku>( r ) oc 
exp(iur/c) describes a wave with frequency \u\ carrying constant energy flux. Within 
the cavity, the function satisfies the Helmholtz equation with the complex wavenumber 
Qk{w)/c. For a given u, there exists an infinite set of eigenmodes labeled by the discrete 
(multi-)index k. The functions V , fc^( r ) are the right eigenfunctions of a certain non- 
Hermitian operator. The conjugate left eigenfunctions, (f>k u (r), behave as exp(— iuir/c) 
in the asymptotic region and have eigenfrequencies [f2 fc (u;)]* inside the cavity. The CF 
eigenfunctions ipku>( r ) and 4>ku}{r) form complete biorthogonal bases of the open cavity. 

Let us represent the electric field inside the cavity as a superposition of 
eigenfunctions ipku{ r )- After substitution of the expansion (A. 6) in the left-hand side of 
Eq. (18), application of differential equation (A. 2) and biorthogonality relation (A. 4), 
we obtain an equation for the amplitude of the mode k, 

[Q 2 k (uj) - uj 2 ] a ku] = 

47T ^ 1 f 

where 0fc qw is the Fourier transform of 0fe w (r), defined as in Eq. (9). The sum Y^=i can 
be substituted by the average Af(- ■ • ) over the positions r° and velocities Vj of active 
atoms. The distributions of the positions and velocities are uncorrelated: the former are 
distributed uniformly within the cavity and the latter obey the Maxwell distribution. 
We approximate 

(e iq ' r ?) „ * \ I dre^ = ^fe(q), (20) 

where (27r) 2 [5j(q)]* is the Fourier transform of the characteristic function of the cavity 
interior and A is the cavity area. <5j(q) has a 5-function property f <iq5j(q) / q = / q= o, 
if / q is the Fourier transform of a function /(r) that vanishes outside of the cavity. By 
substituting the linear and third-order contributions to the dipole moment, Eqs. (14) 
and (16), we transform Eq. (19) to 

H^-^a^^V^+vit (21) 
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where the linear and third-order terms are, respectively, 

x ((u- q -v) 2 D(u-q-v)) , (22) 



V ^ =ic(3) J^8 J dqdqd^'du'duj"uj'D(uj r ) 0£ qa ,£ q _ q ,_ q / V / 

x ^(u - q • v) 2 D(u - q • v) D\\[u - u" - (q' + q") • v] £7 q # w _ w ,_ w ##_ q //. v £ q / V+q ». v ^ 

(23) 

In the velocity averages we dropped the index j labeling individual atoms. The 
coefficients are given by 

m t AN Af d 2 

c (1) = 4tt " - , 24 
^ AN J\f I f d 2 \ 2 

c ( 3 ) = 8tt —- . (25) 



y/eA hv)\\ \hq 

The field amplitudes in Eqs. (22) and (23) can be expanded in the CF quasimodes 
as in Eq. (A. 6). Therefore, Eq. (21) represents a system of algebraic equations for the 
quasimode amplitudes a^ w . 

2.3. Lasing modes 

The frequency integrals in Eq. (23) can be eliminated by approximating the electric field 
in the time domain with a superposition of constant-frequency oscillations, the lasing 
modes. It is well known that when the pumping increases above a certain threshold 
value, the first lasing mode appears. For higher pumping levels, additional modes start 
oscillating. In the multimode regime, the procedure of representing the field spectrum 
as a finite number of frequencies is an approximation: to preserve self-consistency of the 
laser equations, the oscillations at beat frequencies have to be neglected. 
We start with an ansatz 

E^ = 2nJ2 [ E H S(u ~ wj) + 5(u + w,)] , (26) 
i=i 

which assumes the existence of N m lasing modes with frequencies uji > and guarantees 
the realness of E(r,t). By expanding £j q in the CF quasimodes as 

Eiq = e~ 1/2 ^ aik^kqux , (27) 

k 

we obtain the lasing- mode representation for the coefficients a^, 

a k u; = 2tt^2 \ a ik S(u- ui) + a* lk 8(u + wj)] . (28) 
i=i 
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Usually, the quasimode with the longest lifetime starts oscillating at the lasing 
threshold. At higher pumping levels, the nonlinear terms in the laser equations mix the 
quasimodes, and the lasing mode no longer coincides with a single quasimode. Thus, 
the expansion (27) is nontrivial, as it contains, in general, more than one term. 

After substituting the mode representations (26) and (28) in the laser equations (21) 
and collecting the terms proportional to 5(uj — u>i) [or, equivalently, 5(ou + u)i)], we would 
obtain a rather cumbersome equation for the coefficients a^. We write this equation 
in a simplified form by applying the standard approximation 5u <C is, where 5u is 
the frequency spacing between the lasing modes. This approximation corresponds to 
neglecting the second time derivatives in Eqs. (1) and (2). Additionally, we neglect the 
terms containing D\\{2v) ~ 7||/V <C 5u)/v. The resulting equation reads 

(n k -Ui)a lk = PS t ) +Vg ) , (29) 

with 

T>$ = -^(2^ / dq^Wwi-q-v)),^, (30) 
V% = i C <V£^ J W^W,, 

x {(D(ui - q • v) [D(u v - q" • v) + D*(u v - q' • v)] L>||[(q' - q") • v]) v £, >q+q ,_ q „£,, q „ 
+ (1 - 5 w ){D{ui - q • v) [D(ui - q" • v) + D*(u v - q' • v)] 

x D\\[wi - u v + (q - q") • v])v^i' iq +q'- q »Siq»} • (31) 
In the equations above we introduced the notation fl k = fifc(z/), fcq = (j)^, and 

%)=(l-»Y) ( 32 ) 

which approximates D(u) for positive frequencies. The field amplitudes are assumed 
to be expanded in the CF wavefunctions ^ q = V'feq^j as i n Eq. (27). The nonlinear laser 
equations (29) can be solved for the frequencies and amplitudes of the lasing modes, 
e.g., by using the iteration procedure described in Ref. [19]. 

The term proportional to (1 — 5w) in Eq. (31) is usually omitted in the case of fixed 
atoms, v = 0. In such a system, the factor D\\{uji — u^) ~ ~ 7||/t± *C 1, arising 

due to oscillations of the population inversion, makes this term small compared to the 
first term in the braces. In Sec. 3.2 we argue that for a gas laser the second term is 
negligible as well. 

Although the algebraic equation (29) has been derived above for the steady-state 
lasing oscillations, exp(— iuit), it can be easily generalized [14] to the case of slowly- 
varying amplitudes a^(i) by adding the term — ia\ k to the left-hand side. The resulting 
first-order nonlinear differential equation can be used to study stability of the lasing 
modes. 
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3. Properties of the laser equations 

In this section, we simplify the right-hand side of the laser equations (29) by taking into 
account the resonant dependence of the quasimode wavefunctions on q. In particular, we 
show that the linear gain-induced coupling of quasimodes (Sec. 3.1) and the nonlinear 
contribution proportional to (1 — 8u>) (Sec. 3.2) are negligible. 



3.1. Linear coupling of quasimodes 

After the field amplitude is expanded in the basis functions the linear 

contribution (30) can be written as 



V 0) _ 

' Ik ~ 



The matrix elements 



V kk '(ui) = — ^ j dq0fe q (D(wj-q-v)) v ^ fc / q 



(33) 



(34) 



describe possible linear coupling of the quasimodes via the Doppler shift in the frequency 
dependence of the linear gain. A similar coupling in the coordinate space arises in 
the systems with inhomogeneous refractive index or pump power distribution [11, 
14]. Below we argue, however, that the off-diagonal elements Vkk'i^i) are negligible, 
i.e., the quasimode wavefunctions correctly describe the lasing modes in the linear 
approximation. 

First, we evaluate the velocity average in Eq. (34). In the Doppler limit 



\u>i — v\ ~ uv/c ^> 7j_, 
typical for the gas lasers, the function (32) can be approximated as 

1 



D{uj) ~ 7j_ 



tt5(Acj) -HP 



Auj 



Acu = uj — v, 



(35) 



(36) 



where P denotes the principal value. Then the angular average is given by 

(D(u-qv cosy?))^ 



7± 



(qv-\Au\) .6(\Au\-qv) 



a/ q 2 v 2 — Auj 



a/ Auj 2 — q 2 v 2 



(37) 



where 9(x) is the Heaviside step function. After the averaging over the two-dimensional 
Maxwell distribution, 



W(v) alv 



^e- v2/v2 vdv, 



we find 



(D(uj - q • v)>, 



qv 



TTV 



_J_ e -(Au>/qvf + *_ F 
2y^7T 7T 



AuA 
qv J 



(38) 



(39) 
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where F(x) = exp(— x 2 ) J* exp(y 2 ) dy is the Dawson integral [20]. As a function of q, 
(D(u — q • v)) changes on the scale Aq ~ |Au;|/U ~ v/c, for typical Au. 

Next, we note that ip kci and <p kn are significant only when the magnitude of q 
lies within the band \q — q k \ % Sq k , where q k = \feReQ k /c and 5q k = — y/eImQ k /c. 
Therefore, the matrix elements V kk >(oJi) vanish if \qk — q k >\ ^> Sq k . In the reverse case, 
\<lk — Qk'\ ^ Sq k , w e can pull (D(u>i — q • v)) v out of the integral (34), since this function 
is constant on the scale Sq k <C Aq. Then Eq. (A. 5) yields the diagonality property, 

V w (ui) = (Dfa - q • v)) v?=go 5 kk ,, (40) 

where the substitution of q = qo = \ftv /c is consistent with the level of approximation 
in Eq. (29). We note that the linear-gain frequency dependence, given by the real part 
of Eq. (39), has a Gaussian shape, similar to one-dimensional gas lasers [see Eq. (B.3)]. 

3.2. Nonlinear terms 

The resonant dependence of the wavefunctions tp kci and 0fc q on q, mentioned previously, 
makes it possible to simplify the nonlinear terms of Eq. (29) as well. The condition 
\q ~ Qk\ <5<?fc implies one of the following restrictions on the integration domain in 
Eq. (31): (i) |q' - q"| < 5q k , (ii) |q - q"| < 5q k , or (iii) |q + q'| < 5q k . Only 
in the first case we can estimate -D||[(q' — q") • v] ~ min(l, / -f\\/5q k v). In the cases 
(ii) and (iii), we have -D||[(q' — q") ■ v] ~ r y\\c/uv <C mm(l,7\\/5q k v). Similarly, 
D\\[u)i — OJy + (q' — q") • v] <C min(l, ^w/Sq^v) in all three cases. Thus, we can neglect 
the terms in Eq. (31) that are proportional to (1 — 8u>) and rewrite this equation with 
the new integration variables Q = (q' + q")/2 and 5q = q' — q" corresponding to the 
choice (i): 

x D(u i — qov cos (p) C(u i> — qov cos (p') ) , (41) 

/ v 

where C{oj) = ReD(u>), and q = Q = q is set in all nonresonant functions. The value 
of the correlator 

2 roo poo p 

y( ¥ ,,^,v)=^ I jf dqj dQ J ^q)L>||(5q.v) 

X 0fcq^,q+5q £, r',Q+<5q/2^',Q-5q/2 (42) 

is system specific. Here (p (<//) is the azimuthal angle of the vector q (Q) counted from 
the direction of v. We can approximate D\\(5q • v) ~ 1 under the assumption of long 
living quasimodes, such that 5q k <C "f\\/v. In this case, the correlator depends only on 
the direction of v, but not on its magnitude. The absence of functions D\\ from the laser 
equations is an indication that the approximation of constant population inversion is 
employed. 



Gas lasers with wave-chaotic resonators 



11 



4. Chaotic gas laser 

4-1. Nonlinear terms in the laser equations 

We apply the general theory to a gas laser with wave-chaotic resonator. The notion of 
a chaotic system implies that the shape of the resonator does not permit a separation 
of variables in the wave equation and the relevant modes have a typical wavelength A = 
27r/go "C L. Famous examples of chaotic resonators are the Sinai and stadium billiards 
and their desymmetrized versions [21]. The classical ray dynamics, defined in the short- 
wavelength limit, is assumed to be ergodic, i.e., the generic ray trajectories fill densely 
the available phase space. This assumption sets a limitation on the degree of openness 
of the resonator: the lifetimes of the CF modes contributing to the lasing modes are to 
be longer than the time needed for the ray trajectories to establish ergodicity. This time 
is of the order of inverse Lyapunov exponent. We will consider a weakly open resonator, 
such that its eigenmodes are well approximated by real eigenf unctions ipk{ r ) ~ 0fc( r ) 
independent of the parameter 00. 

Let us examine a product of two wavefunctions evaluated at nearby points in the 
wavevector space: 



Since the relevant magnitude of 5q is much smaller than the optical wavenumber go, 
the function e l<5q ' R changes on the length scale much longer than A. On the other hand, 
the wavefunction product self-averages in the chaotic cavity over a correlation area of 
order A 2 [15]. Consequently, the product can be substituted with its average value, 
expressed in terms of the correlation function 



where the wavefunctions are assumed to be normalized. Different eigenmodes, k 7^ k', 
are uncorrelated in the chaotic system, whereas the autocorrelation function in two 
dimensions is approximated by a Bessel function [15], 



This expression shows that the statistical characteristics of the eigenstates are isotropic 
(this approximation fails within a wavelength of the cavity boundary). Scars [22], i.e., 
the excess probability density around unstable short periodic ray trajectories, should not 
significantly influence the correlation function. Indeed, the excess probability density 
depends on the inverse Lyapunov exponent, which scales with the system size. On 
the other hand, the fraction of the phase space occupied by a scar reduces with the 




(43) 




R 



(44) 



Ckk(Sr) = J (q k 5r). 



(45) 
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wavelength, which, in gas lasers, is several orders of magnitude smaller than the system 
size. The product (43) is proportional to the Fourier transform of the autocorrelation 
function, (2n/qk) 5(Q — q k ), which makes 

^fc,Q+<5q/2 ^fc',Q-5q/2 = ^j- 5(Q - q k ) 5 x (5q). (46) 

We use this result to calculate the third-order contribution (41) to the laser 
equations. An important simplification that comes with the stochasticity of the 
wavefunctions is the isotropy of the correlator (42). Hence, this function can be pulled 
out of the angular integrals, and we find 



v=i \ k 



( 3 ) = ^1,(3)„,, \^(\^| a/ , fc/ |2 



x ((D(ui- q v cos (p)) (C(u v -q v cos (p)) ) . (47) 



' v 



4-2. Single-mode regime 

We apply the laser equations (29), with the right-hand side given by Eqs. (33), (40), 
and (47), to determine the dependence of the mode intensity on the frequency in the 
single-mode regime. If the pumping level is not far above the threshold, only one 
quasimode will contribute to the lasing mode. Consequently, a single term with I' = I 
and k' — k remains in the summations in Eq. (47). Taking the imaginary part of 
Eq. (29), we obtain an expression for the intensity, I k = \ai k \ 2 /e, 



(C(v-q v cosp)) 2 



(^(C(uj — qov cos (p)} 

where the intensity unit 



' V 



.(a^\ 2 AAV 

\Q0 V J — 

AA 



(48) 



7 °=8^?^(^-^ cos <)r ^ 

is defined by the condition h(v) — > h, as AA — > oo. The threshold value of the 
pumping strength AA at u> = v is given by the parameter 

a at eAhq v 

ANo > k = Kk ^M^ (50) 

where n k = — Im Q k is the decay rate of the quasimode. 

When evaluating <^(£)^ , it is not possible to use the approximation (37) for the 
angular average, because the velocity integral would diverge. The angular integral can 
be calculated analytically for a finite value of •y±/qoV, however, the resulting expression 
is rather cumbersome. The remaining integral over the distribution W(y) has to be 
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Aw/q Q v(2D), A<j/V2~ q Q v(\D) 

Figure 1. Relative intensity of the single lasing mode for the two-dimensional (2D) 
chaotic laser [Eq. (48), thick curves] and a one-dimensional (ID) laser [Eq. (B.7), 
thin curves] as a function of the frequency separation Aw from the atomic-line 
center. The scaling of Aw for these systems differs by a factor of \[2 (see the text). 
Parameters: relative pumping strength ANq/ AiVo.fc = 1-3; relative homogeneous 
linewidth r yj_/qo'v = 0.01 (solid curves), 0.04 (dashed curves), and 0.07 (dash-dotted 
curves). 



evaluated numerically. In the special case u — u the result can be expressed in terms of 
a special function: 

(C{v - q v cos = |- e x E 1 {x) ) (51) 

\ 2 

71 



x = , (52) 

\QovJ 

where Ei(x) = f^°dy e~ y /y [20]. 

The relative intensity Ik(co)/Io is plotted in Fig. 1 for several values of j^/qov in 
the Doppler limit (35). The frequency dependence is characterized by a local minimum 
of intensity at the center of the gain curve, u) — v. The frequency interval where the 
intensity is reduced has a width of the order of go^- Formally, this follows from the 
analysis of the denominator of Eq. (48). In fact, the integral 

( C(lu — q v cos (p)) in = — [ — xo . 9 (53) 

X V ,llp ttJo l + (Au;-go^cos^) 2 / 7 i v ; 

receives its major contribution from the angles in the interval 5(f ~ 7±/<?o^ around 
ip = cos -1 (Aoj/q v). Therefore, the integral changes its order of magnitude from , j±/qoV, 
for |Ac<;| < q v — 7_i_, to (7j_/Ao;) 2 , for |Aa;| > q$v + -y±, i.e., go^ is the characteristic 
frequency scale. 

The results for the chaotic laser are compared in Fig. 1 with the relative mode 
intensity in a one- dimensional gas laser (Appendix B). In the former (latter) case, Au is 
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normalized by q$v (y/2qov). Apart from the visual convenience, the difference in scaling 
is motivated by the relation, v 2 d = v^2^id, between the root mean square velocities 
in two and one dimensions at a given temperature. Lasers with uniaxial resonators 
exhibit a much narrower intensity minimum, known as Lamb dip [1]. Its width is of 
the order of the homogeneous linewidth 7_|_ <C qov. In order to understand the origin of 
the dip, it is convenient to treat the lasing mode as a superposition of left- and right- 
moving waves. The atoms moving with the velocity v x ~ Au/q x effectively interact 
with only one of the two waves, depending on the sign of Au and v x . However, when 
\Au\ < 7j_, the atoms at rest interact with both running waves. This leads to a stronger 
gain saturation near Au = than away from the line center. In the two-dimensional 
chaotic laser this effect is washed out, since a mode with the frequency oj draws its 
gain from a larger ensemble of atoms with velocities v > \Au\/qo. In particular, let us 
model the chaotic wavefunction locally by a random superposition of partial plane waves 
with wavevectors of fixed magnitude y/eu/c and isotropically distributed directions [16]. 
Then, all partial waves of the central mode u — v shall interact with the atoms at rest, 
but each partial wave interacts, in addition, with atoms moving perpendicular to its 
direction of propagation. Thus, the minimum of intensity results from a more subtle 
interplay of linear amplification and nonlinear saturation on the same frequency scale. 

The results obtained above for a wave-chaotic resonator cannot be directly applied 
to resonators with mixed or regular classical ray dynamics, where the hypothesis 
of random, isotropically distributed, and uncorrelated eigenstates is no longer valid. 
Nevertheless, one can conjecture that increasing dimensionality of the system will make 
the intensity minimum shallower and wider (its depth and width being defined relatively 
to the intensity maximum and the inhomogeneous linewidth, respectively), but will 
not lead to its complete disappearance. As a practical consequence, the frequency 
stabilization [4] and spectroscopy [5] might still be possible in imperfect quasi-one- 
dimensional resonators, where the lasing mode deviates from a superposition of two 
counter-propagating plane waves. It is worth mentioning that in both uniaxial and 
chaotic resonators the relative depth of the minimum decreases with increasing 7^. 

5. Conclusions 

In the present article, the semiclassical laser theory was extended to gas lasers with 
two-dimensional resonators of arbitrary shape. It was shown that the linear coupling 
between the modes, resulting from the Doppler shift in the gain frequency dependence, is 
negligible. The nonlinear coupling was considered in the third order of the perturbation 
theory in the electric field. The nonlinear terms in the laser equations that arise due 
to pulsations of the population inversion were identified and neglected. The criterion 
(v/c)v ^> 7||, allowing one to approximate the population inversion as constant, is 
usually well fulfilled (is is the atomic-transition frequency and v is the average velocity). 

The general theory was applied to lasers with two-dimensional weakly open 
resonators of wave-chaotic geometry. In the short- wavelength limit, the eigenf unctions 
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of the resonator can be locally approximated by a superposition of plane waves with a 
fixed wavelength propagating in random directions. The applicability condition for this 
description can be estimated as c/y <C i« € I, where L is the size of the resonator 
and w is the width of the output window. The isotropic statistical characteristics of the 
eigenfunctions lead to decoupling of the field distribution from the Doppler-broadened 
gain curve in the nonlinear terms of the laser equations. The single-mode intensity, as 
a function of the laser frequency, has a local minimum at the frequency of the atomic 
transition. The width of the minimum scales approximately with the inhomogeneous 
broadening (v/c)v and only weakly depends on the homogeneous linewidth j±. This 
property distinguishes the intensity minimum in chaotic resonators from the Lamb dip 
in uniaxial resonators, which has the width of 7j_ <C (v/c)v. 
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Appendix A. Constant-flux quasimodes 

We consider an open cavity with a real dielectric constant e(r) = const inside and e = 1 
outside of the system's boundary. The constant-flux (CF) modes [12] ip k u){ r ), depending 
on the parameter oj > 0, are defined as solutions of a non-Hermitian eigenvalue problem. 
Explicitly, the functions must satisfy the differential equation 



in the exterior of the cavity with the outgoing-wave boundary conditions at infinity, 
while inside the system the same mode satisfies a different equation: 



For any given value of u, the complex eigenfrequency fl k (ou) is quantized [k is a discrete 
(multi-)index labeling the modes], because the solutions are required to match smoothly 
at the interface. 

The conjugate wavefunctions (p kL0 (r) obey Eq. (A.l) with the incoming-wave 
boundary conditions outside the cavity and the equation 



inside the system. The CF quasimodes and their conjugates are biorthogonal, and can 
be chosen to satisfy the condition 



-c 2 v 2 i) kul = u 2 i) kul 




(A.2) 




(A.3) 




(A.4) 
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where the integration is over the interior I. A similar relation, 

J^y f d< l^fcq«^fc'qw = 5 k k>, (A.5) 

is valid in the q representation, defined as in Eq. (9). Additionally, the wavefunctions 
possess the properties: fcw (r) = [^(r)]* and fcqw = ^ _ q>w . 

A Fourier component -E^r) of the lasing field can be expanded in the CF modes as 

E u (r) =e- 1/2 ^a fcw ^(r), (A.6) 
k 

a^^e 1 ' 2 J drPtoWE^r). (A.7) 

When continued to the exterior, this expansion yields a wave at the frequency u 
propagating in the free space away from the system. In the stationary lasing regime u 
becomes restricted to a finite number of values corresponding to the frequencies of lasing 
modes. To guarantee the realness of the field E(r,t), we need to extend the definition 
of the CF quasimodes to negative frequencies by 

^fc,-«(r) = 0*w(r), <£k,-w(r) = ^(r). (A.8) 

Since the eigenfunctions <t> ku) {r) depend only on the squared complex frequency VL 2 k {oj), 
the frequencies VL k {oj) are determined up to an overall sign. To describe the field 
distribution with the outgoing flux and no incoming flux, the sign is chosen according 
to the requirements: 

Re[Q k (u)} > 0, lm[Q k (u)] < (u > 0), (A.9) 
n k (-w) = -p k (u)]\ (A.10) 

In connection with the definition (A.8), this choice provides for negative Im[fi fc (u;)] for 
any u>, as is expected from the decaying quasimodes. 

Appendix B. One-dimensional gas laser 

We apply the general theory to a one-dimensional gas laser that was studied in detail in 
Ref. [1]. Specifically, we consider a resonator defined by e = 1 in the region < x < L, 
with a perfect mirror at x = and a weakly transmitting mirror at x — L. The 
quasimode eigenfunctions have the form 



piClk-x/c —iQf,x/c 

M*) = -T^r ~ —jwt = 4 + \x) + 4-\x), (b.i) 



where we approximate q k = Refl k /c ~ nk/L (k = 1,2,...) and assume 5q k = 
— lmQ k /c C In principle, Sq k can be determined by modeling the semitransparent 
mirror as a dielectric layer with large dielectric constant and solving Eqs. (A.l) 
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and (A. 2), but the specific value of 5q k is not important for our purposes. In the q 
representation, the functions 



2L qk-iSq k Tq 



(B.2) 



have a resonant dependence on q, peaked at ±q k and having a width of Sq k . 

As in the two-dimensional case, the velocity average of the homogeneous gain curve 
can be computed within the approximation (36). Averaging over the one-dimensional 
Maxwell distribution yields 



1 /Aw 

2 yqV 



(B.3) 



The third-order contribution is obtained from Eq. (31), adapted to the one-dimensional 
case. The electric field can be represented as a superposition of the quasimode 
wavefunctions (B.2). The velocity averages depend on q on a scale larger than Sq k , and 
can be pulled out of the integrals and evaluated at the resonant values of q. The integrals 
over the wavefunctions have a quite simple form in the coordinate representation: 



<Pk",q' 



J dqdq'dq"^l q ^' ) q+ql _ qll 



rk'", q " 



(B.4) 



where a, a', a", a'" = +, — and we took into account that <j>%{x) = ^ k {x). The phase of 
the integrand vanishes when the indices obey the resonance condition 



ah + a'k' - a"k" + a"'k"' = 0. 



(B.5) 



In the single-mode near-threshold regime, when only one quasimode is excited, six of 
the 16 possible combinations of indices are resonant. However, only four combinations, 

(a, a', a", a'") = (H h+), ( — I ), (H ), ( — yield q' = q" at resonance, 

which results in D\\ — \. In the limit 5q k L <^ 1, the resonant value of the integral (B.4) 
is equal to 1/4L. The velocity averages in Eq. (31) with the resonant values of q — ±q k 
and q' = q" ± q k sum up to 4 (D(uu — q k v) [C(uj — q k v) + C(u + q k v)]) v - This quantity 
was calculated in Ref. [1]. The real part in the Doppler limit is given by 



(C(u - qv) [C(u - qv) + C(u + qv)]) v 

2' 



IT 7 ± 
= \ ~ — eX P 

qv 



qv 



[1 + £(«>)] 



(B.6) 



The imaginary part of Eq. (29) yields an expression for the mode intensity I k = \a ik \ 



= 2 



Ay 
2 \ q v 



1 + C(u>) 



(B.7) 
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where we approximated q^ ~ q . The parameters 



/o = -2^T' (R8) 
AiV 0fc = * qkfWL - (B.9) 



are defined in Sec. 4.2. 
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